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A general localization mechanism for waves in dissipative 
systems is identified that does not require the bistability of 
the basic state and the nonlinear plane- wave state. The mech- 
anism explains the two-dimensional localized wave structures 
('worms') that recently have been observed in experiments on 
electroconvection in nematic liquid crystals where the transi- 
tion to extended waves is supercritical. The mechanism ac- 
counts for the propagation direction of the worms and certain 
aspects of their interaction. The dynamics of the localized 
waves can be steady or irregular. 
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A striking feature observed in a number of pattern- 
forming systems with large aspect ratio is the sponta- 
neous localization or confinement of the pattern to a 
small part of the system although the system is trans- 
lationally invariant. Presumably the best studied struc- 
tures of that type are the quasi-one-dimensional pulses 
of traveling waves that have been found in convection in 
binary mixtures (e.g. [jD). They have been described the- 
oretically as perturbed solitons (e.g. ||||) and as bound 
pairs of fronts . Other quasi-one-dimensional local- 
ized structures have been found in Taylor vortex flow 
, directional solidification g , cellular flames || , and 
in models for parametrically driven waves JXG ] ■ 

In two dimensions localization appears to be harder 
to obtain; investigations of binary mixture convection 
have only led to long-lived, but eventually unstable 
wave pulses jnj. Only recently have truly stable two- 
dimensional localized structures been found in paramet- 
rically driven surface waves in granular media and in a 
highly viscous fluid |Q. In most systems the localized 
structures arise in a regime of bistability. In such situa- 
tions they can often be considered as a pair of fronts that 
are bound to each other by dispersion f4j , an additional 
mode H , or non-adiabatic effects |l3|-|16[]. 

Very recently [ p"7| a two-dimensional localized wave- 
state has been found in a system in which the transi- 
tion to the extended waves is supercritical, i.e. in which 
there exist no fronts connecting the basic and the nonlin- 
ear state. This rules out the mechanisms of localization 
mentioned above. Thus, previous efforts to explain these 
states turn out to be insufficient |l5| and the origin of lo- 
calization has remained quite puzzling. The experiments 
have been performed in electroconvection of nematic liq- 
uid crystals and due to their appearance the new local- 
ized states have been called 'worms'. 



In this paper we discuss a general mechanism that can 
lead to localization even if the bifurcation to the extended 
waves is supercritical. It is based on the presence of an 
additional, weakly damped field that is advected by the 
waves and that in turn affects their growth rate. The 
mechanism explains the shape of the worms, their direc- 
tion of propagation and certain aspects of their interac- 
tion. 

Before the theoretical model is introduced a number 
of features of the experimentally observed worms need to 
be discussed. Due to the preferred direction associated 
with the liquid crystal the system has axial anisotropy. 
In the regime in question convection arises in the form 
of waves that travel at a fixed angle relative to the axis 
of anisotropy. Due to reflection symmetry there are four 
such directions of propagation (left- and right-zig, and 
left- and right-zag) as indicated in fig.[j]. The worms con- 
sist either of a combination of the left-traveling waves or 
of the right-traveling waves. 
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FIG. 1. Sketch of the orientation of the group velocities 
ua,b and of the directions of advection Iia.b (cf. (|l|- ||) be- 
low). 



According to recent calcuations the transition to ex- 
tended waves is supercritical in this regime This is 
consistent with the experiments, where for larger con- 
ductivity of the liquid crystal extended chaotic waves 
are observed to arise supercritically and no indication 
of a change to subcritical behavior of these waves is seen 
|]l7f . Nevertheless, the worms arise already at parame- 
ter values at which the state without convection is stable 
Jl9| . In the x-direction the worms are quite long and 
their length appears to be somewhat variable. In the y- 
direction, however, they are very narrow and their width 
is fixed. Unless perturbed, the worms travel only in the 
x-direction. In that direction the envelope of the worm 
has a characteristic shape: it rises rapidly to its maxi- 
mum at one end and decays slowly over the length of the 
worm. The worm travels toward the end with the large 
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amplitude (its 'head'). 

The worms occur very close to threshold and the bi- 
furcation to the extended traveling waves is supercriti- 
cal flgl . Therefore a set of coupled complex Ginzburg- 
Landau equations is considered. Since in a given devel- 
oped worm a single set of zig- and zag- waves is observed 
jl7j we consider only these two waves. 

Within the Ginzburg-Landau equations no sustained 
waves are possible below threshold since all nonlinear 
terms are damping. We extend therefore these equations 
introduce an additional, weakly damped mode. The mo- 
tivation for this mode arises from previous work in the 
context of binary-mixture convection where it was found 
by one of the authors that pulses of traveling waves can 
arise even for a supercritical bifurcation if the waves ad- 
vect a mode that feeds back into their growth rate (2^] . In 
addition, the weak-electrolyte model, which agrees quan- 
titatively with the experiments with respect to the linear 
properties suggests that a charge-carrier mode be- 
comes slow in the regime in which worms appear 18 1. 

A minimal model for the advection of a scalar mode 
by zig- and zag- waves is given by 



b x d 2 A 



b y d 2 y A 



2ad 2 xy A (1) 



d t A= -UA-VA + llA 1 '■ >' 

+fCA-c\A\ 2 A-g\B\ 2 A, 
d t B = -vlb -VB + LiB + b x d 2 x B + b y d 2 y B - 2ad 2 y B (2) 
+fCB -c\B\ 2 B - g\A\ 2 B, 



d t C = 5d 2 x C -aC + h A - V\A\ 2 + h B ■ V\B\ 



(3) 



The equations for the complex wave amplitudes A and 
B are the usual complex Ginzburg-Landau equations for 
oblique waves |2^ ]. The equation for the scalar mode 
C is obtained by considering the currents = h^yl) 2 
and ]b = ^-b\B\ 2 that are due to the advection by the 
respective waves. In addition, damping and diffusion of 
the mode C is allowed. In the following we first focus on 
the effect of C and will mostly consider the dispersionless 
case in which all coefficients are real. 

Eqs.(|]-^) are solved numerically using a pseudospec- 
tral code with an integrating- factor/Runge-Kutta time- 
stepping scheme. Fig.^| shows a worm solution obtained 
in this way. For illustration purposes the top part shows 
the real part of the quantity Ae l<1A ' T +Be l<lB ' r which gives 
an indication of how this solution would appear in exper- 
iments. The mode C is shown in the bottom part. As in 
the experiment the convective amplitude is large at one 
end of the worm ('head') and decays slowly towards the 
other end. In agreement with the experiments, the worm 
propagates toward their head. 

To understand how this localized solution can arise al- 
ready below the threshold /i = although the transition 
to extended waves is supercritical, two one-dimensional 
reductions of are considered. The long and narrow 
shape of the worms suggests that the localization mecha- 
nisms along and transverse to the worm differ from each 
other. If one ignores the ^-dependence, (|l|-|^) reduce to 
two equations describing standing waves coupled to C 



and the worm is replaced by a stationary pulse of stand- 
ing waves. Such a solution is shown in fig.^. It exists 
because the C-mode enhances the growth rate in the cen- 
ter of the pulse and keeps the two components A and 
B together. Since the currents generating C vanish if 
|A| = \B\, the pulse can only exist if the two traveling- 
wave components are shifted with respect to each other. 
The shift is due to the group velocity in the y-direction. 
Thus, the worm disappears in a saddle-node bifurcation 
when the group velocity is reduced below a certain value. 
The standing- wave pulse is stable since the production of 
C increases with the shift between A and B. 





FIG. 2. Numerically determined worm-solution of (y-|3|) for 
u A = (0.05,0.75), fj, = -0.025, b x = b v = 1, a = 0, c = 1, 
g = 0.5, / = 1, a = -0.02, S = 0.7, h A = (1, 1). a) gray-scale 
plot of Ae i<1A ' T + Be iqB ' r . b) mode C. Note the depression 
of C along the sides and the back of the worm. The worm 
travels to the left (toward its 'head'). 
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FIG. 3. Standing-wave pulse in one dimen- 
sion (y-direction). Parameter values corresponding to flg.g. 

For the understanding of the stability of the worm so- 
lution in fig.§ it is crucial to note that the standing-wave 
pulse shown in figj^ arises in a subcritical bifurcation al- 
though the bifurcation to the extended standing (or trav- 
eling) waves is supercritical. Thus, the standing-wave 
pulse can exist already for /i < and can coexist with 
the basic state A = B = C = 0. 

The standing-wave pulse can be considered as a build- 
ing block for the worm: in terms of a one-dimensional 
reduction along the x-direction the worm can be seen as 
a pair of fronts connecting the basic state with the coex- 
isting standing-wave-pulse state. For stability these two 
fronts have to exhibit a repulsive interaction. With re- 
spect to the x-directon the zig and the zag making up the 
worms are waves traveling in the same direction. Based 
on results for traveling- wave pulses QJ^,^] it is therefore 
expected that the fronts can interact ij' via the wavenum- 
ber which is driven by the dispersive terms, and ii) 
via the additional mode C ]5|,p3|]. 

If the interaction via C dominates, a simple connection 
between the stability of the traveling-wave pulse (here, 
the worm) and its direction of propagation emerges || . 
As shown in figj^b there is a large peak of C at the head 
and a shallow depression of C at the tail. Since positive 
(negative) C enhances (reduces) the local growth rate of 
the convective mode both fronts are pushed to the left by 
C. If the worm travels towards its head, the depression 
at the trailing front will be reduced by the remnants of 
the positive peak ahead of it. This is not the case for the 
leading peak. Thus, the trailing front is pushed less to 
the left than the leading front amounting to a repulsive 
interaction. For the traveling- wave pulses this interaction 
has been determined previously in the limit of weak dif- 
fusion 8 and damping a and for widely separated fronts 

For the present case of fronts in the standing-wave 
amplitude the corresponding calculation would be con- 
siderably more involved. We expect, however, that the 
same qualitative picture holds. This would imply that 
the worm is stable if it travels towards its head and un- 
stable otherwise [||. Indeed the experimentally observed 
worms travel towards their head. It should be noted that 
sufficiently strong dispersion can lead to additional sta- 
bilization of the worms |2^] . 

When increasing [i < the width of the worm and 
the maximal value of C remain essentially unchanged, 
but the worm grows in length. Consequently, the spa- 
tial integral Af = J \A\ 2 dxdy, which corresponds to a 
kind of Nusselt number, increases smoothly except for a 
very small jump in Af at the saddle-node bifurcation in 
which the worm first appears. In a sufficiently large sys- 
tem, when the threshold fj, = is surpassed a sequence 
of transitions to more than a single worm occurs. When 
the worms become long enough to span the whole sys- 



tem they loose their pronounced head structure and C 
becomes independent of x. 

In the experiments the worms exhibit a typical spacing 
in the y-direction which decreases with increasing applied 
voltage. To investigate this aspect within (Q-^), the tem- 
poral evolution starting from random initial conditions 
is followed. Fig.^a gives the y-position of the emerging 
worms as a function of time for /i = 0.05. Fig.^jb shows 
the C-field of the worms at the final time of the run. As 
in the experiment, the worms approach roughly a typical 
distance in the y-direction and almost travel on 'tracks'. 
When /i is increased the number of 'tracks' is found to 
increase. Fig.^b also demonstrates that the localization 
mechanism in the y-direction is much stronger than that 
in the x-direction: while all worms have essentially the 
same width they vary substantially in their length. This 
is expected since the localization in the x-direction is 
achieved via the interaction of fronts, which is very weak 
and decays exponentially in space. Thus, the timescales 
over which the worms reach their final length is much 
longer and perturbations from other worms affect the 
length more strongly than the width. 




Time 

FIG. 4. Formation of worms from random initial condi- 
tions. Parameters as in fig.^| except for /i = 0.05 and 
ua = (0.75,0.75). a) y-location of the maxima in C as a 
function of time, b) C-field for t = 800. 

So far the waves have been taken dispersionless. In the 
experiments it is found that with increasing conductivity 
the worm regime crosses over to a regime characterized by 
bursts in the convective amplitude, i.e. patches or blobs 
of large amplitude appear and disappear in an irregular 
fashion, and finally a regime of extended convection ex- 
hibiting spatio-temporal chaos of patches of zigs, zags, 
and rectangles is reached |f7| . This suggests that disper- 
sion becomes increasingly important along this path in 
parameter space and indeed the imaginary parts of the 
coefficients in the Ginzburg-Landau equations for A and 
B (with C eliminated) increase in this direction [|l8| . 

Fig.|5| shows the result of a run with intermediate val- 
ues of the imaginary coefficients. The values have been 
chosen with some guidance from [ p"8| . In this regime the 
worms turn out to be unstable and start to travel not 
only in the x-direction but also in the y-direction. The 
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motion in the y-direction is driven by imbalances between 
the amplitudes of the zig- and the zag-component which 
lead to enhanced advection of C towards one or the other 
lateral sides of the worm. In addition, the amplitude of 
the worms grows and decays as indicated by the size of 
the symbols marking the y-position in fig.||. These strong 
variations in the amplitude may be related to the burst- 
ing seen in the experiments. 
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FIG. 5. Unsteady motion of worms. {/-location as 
a function of time for ua = (0.05,0.05), = 0.05, 
b x = by = 1 - 0.35i, a = 0, c = 1 + 0.5i, g = 0.5 + 0.9i, 
/ = 1, a = -0.05, 5 = 0.4, h A = (1,1). The size of the 
symbols indicates the height of the maximum of C. 

In conclusion, we have introduced an extension to the 
complex Ginzburg-Landau equations for zig- and zag- 
waves that leads in a simple way to localized, worm-like 
solutions although the bifurcation to extended waves is 
supercritical. The localization mechanisms differ in the 
x- and the y-direction: while the worm is spatially homo- 
clinic in the y-direction, i.e. there is only the basic state 
as a fixed point and the standing-wave pulse represents an 
excursion from it, it is heteroclinic in the ^-direction, i.e. 
it connects the basic state with a nonlinear periodic so- 
lution - the standing- wave pulse - via two fronts. Within 
this framework the worms travel towards their end with 
larger amplitude, in agreement with experimental obser- 
vations. The observed typical spacing between the worms 
in the y-direction is interpreted as the distance over which 
the advected field C suppresses convection. If one as- 
sumes that the damping of the additional mode increases 
with the conductivity the presented model captures qual- 
itatively the experimentally observed change from steady 
worms to extended spatio-temporal chaos via an irregu- 
lar bursting as the conductivity is increased. Extended 
spatio-temporal chaos is obtained when the damping is 
large and the additional mode can be eliminated . 

The localization presented here is similar in spirit to 
that found in parity-breaking bifurcations. There the lo- 
calized structures are drift waves embedded in a station- 
ary, spatially periodic state (rather than an unpatterned 
state) HH|] and the additional field is the wavenumber 
of the underlying pattern ]2j| . 
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